Evidence for isotropic s-wave superconductivity in high-entropy alloys

High-entropy alloys (HEA) form through the random arrangement of five or more chemical elements on a crystalline lattice. Despite the significant amount of resulting compositional disorder, a subset of HEAs enters a superconducting state below critical temperatures, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$T_{\text{c}}<10\,$$\end{document}Tc<10 K. The superconducting properties of the known HEAs seem to suffice a Bardeen–Cooper–Schrieffer (BCS) description, but little is known about their superconducting order parameter and the microscopic role of disorder. We report on magnetic susceptibility measurements on films of the superconducting HEA (TaNb)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$_{1-x}$$\end{document}1-x(ZrHfTi)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$_{x}$$\end{document}x for characterizing the lower and upper critical fields \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$H_{\text{c,1}}(T)$$\end{document}Hc,1(T) and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$H_{\text{c,2}}(T)$$\end{document}Hc,2(T), respectively as a function of temperature T. Our resulting analysis of the Ginzburg–Landau coherence length and penetration depth demonstrates that HEAs of this type are single-band isotropic s-wave superconductors in the dirty limit. Despite a significant difference in the elemental composition between the \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$x=0.35$$\end{document}x=0.35 and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$x=0.71$$\end{document}x=0.71 films, we find that the observed \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$T_{\text{c}}$$\end{document}Tc variations cannot be explained by disorder effects.

www.nature.com/scientificreports/ In this letter, we report on temperature-dependent magnetization measurements of (TaNb) 1−x (HfZrTi) x HEAs at different alloy compositions x to characterize the superconducting state in more detail. Our analysis of the superconducting penetration-depth (T) is in quantitative agreement with BCS theory for an isotropic single-band s-wave superconductor in the weak (to intermediate) coupling limit. Our experimental results further show that, despite the large amount of atomic scale disorder, the observed T c variations for films of different elemental compositions cannot be explained by disorder effects.

Experiment
Films of superconducting (TaNb) 1−x (ZrHfTi) x with nominal x = 0.40 and x = 0.75 have been prepared by magnetron sputtering on the surface of SiN wafers as described in Ref. 13 . In contrast to other HEA types 24 , the zero or small binary mixing enthalpies of the constituent elements in TaNbZrHfTi HEAs favor the formation of a single phase structure, when depositing the HEA film on a substrate held at room temperature 13 . Consistent with our previous study of this compound 13 , the formation of a completely mixed single phase and the absence of other binary phases is confirmed by chemical mapping through energy-dispersive X-ray spectroscopy measurements with the scanning electron microscope. These measurements further facilitate the determination of the actual film compositions x = 0.35 and x = 0.71 that closely match the targeted compositions. X-ray diffraction measurements confirm the single-phase crystallization on a body-centered cubic lattice with a film thickness d ≈ 1 µ m. For details on the binary enthalpies, and the chemical and structural characterization please refer to Sect. A of the Supplementary Materials.
The molar mixing entropy �S = R x i log(x i ) (R-ideal gas constant), of the x = 0.71 alloy, S x=0.71 = − 1.56 R , is comparable to that of the x = 0.35 alloy, S x=0.35 = − 1.53 R . It is interesting to note that the (TaNb) 1−x (ZrHfTi) x alloys with similar mixing ratios x = 0.40 and x = 0.75 were reported to show the highest T c ≤ 7 K and highest H c ≈ 10 T, respectively 13 .We have performed high-resolution DC magnetization measurements using a commercial Quantum Design MPMS3 VSM-SQUID magnetometer under cryogenic conditions T ≥ 1.8 K to characterize the superconducting state of the HEA samples. Measuring their magnetic susceptibility χ(T) , we have determined their T c , the lower H c1 and upper H c2 critical fields as a function of temperature and externally applied magnetic field H.

Results
Zero-field cooling/Field cooling (ZFC/FC) measurement were performed to establish superconductivity in the HEA films, see Fig. 1a,b. Both the x = 0.35 and x = 0.71 film show a diamagnetic response with unity superconducting volume fraction in ZFC measurements at T ≪ T C . The extracted T c = (6.7 ± 0.1) K and T c = (4.3 ± 0.1) K of the x = 0.35 and x = 0.71 film, respectively are in agreement with previous reports 26 . The FC measurements further indicate strong flux pinning. The diamagnetic response of the x = 0.35 film is suppressed by about 60%, whereas the x = 0.71 film exhibits a small paramagnetic Meissner effect 27 .
The penetration depth of a superconductor can be determined through measurements of H c1 and H c2 . We have determined H c1 (T) by mapping out the field response of the HEA films at small external magnetic fields applied in parallel to the film plane. In Fig. 1c, we plot the corresponding superconducting volume fraction − 4π M(H) at different experimental temperatures for the x = 0.35 film (see Sects. B,C of the Supplementary Materials for the corresponding data of the x = 0.71 film and the superconducting volume fraction determination, respectively). At small applied fields, − 4πM(H) exhibits a linear dependence with a slope of unity. This observation is consistent with the diamagnetic response of bulk superconductivity in the HEA films.
The deviation from linearity at larger H occurs at H c1 at which the HEA films enter the mixed phase, i.e., magnetic vortices are penetrating the superconducting volume. We have determined H c1 as the field at which the measured − 4π M(H) data deviate from a linear fit to the small-field region, see Fig. 1c inset. The fitting procedure is described in Sect. D of the Supplementary Materials. The resulting H c1 (T) is displayed in Fig. 1d. While H c1 is strongly suppressed for both alloy compositions at T → T c , H c1 of the x = 0.71 film is about an order of magnitude larger compared to the x = 0.35 film at T ≪ T c . Furthermore, other measurements show that H c1 (T) is not affected by a 45% rotation of H in the sample plane (see Fig. 1d).
We have measured the magnetic susceptibility over a larger field range of − 70 kOe < H < + 70 kOe to further determine the temperature dependence of H c2 . In Fig. 2a, we plot the corresponding − 4πM(H) for representative measurements of the x = 0.35 film (see Sect. D of the Supplementary Materials for the corresponding data of the x = 0.71 film). We observe a significant magnetic hysteresis between forward and backward sweep, indicative of vortex pinning below the irreversibility field H irr (see inset of Fig. 2a). H c2 can be determined from these measurements as the field, at which forward and backward trace deviate from the linear background signal, see marker in Fig. 2a. The resulting H c2 (T) dependence is shown in Fig. 2b. We observe a monotonic, almost linear, decay of H c2 (T) near T c for both alloys.

Discussion
We can accurately describe H c2 (T) by using the Werthamer-Helfland-Hohenberg (WHH) model of conventional superconductors in the presence of spin-paramagnetism and spin-orbit interaction (see Fig. 2b) 25 . Fitting H c2 (T) , we obtain H c2, 0 = (81.8 ± 0.4) kOe and H c2, 0 = (71.9 ± 0.6) kOe for the x = 0.35 and x = 0.71 film, respectively. These values are significantly smaller than the values of the corresponding Pauli paramagnetic limit in the weak coupling limit H P = 18.4T C ( H P in kOe and T C in K) 28,29 . H P = (123.3 ± 1.8) kOe for the x = 0.35 film and H P = (79.1 ± 1.8) kOe for the x = 0.71 film, indicating that superconductivity is rather limited by orbital effects induced by the externally applied field.
We obtain the Ginzburg-Landau (GL) coherence length ξ through the analysis of www.nature.com/scientificreports/ the electron charge. The resulting ξ(T) are shown in the inset of Fig. 2b for both alloy compositions. Their diverging characteristics for T → T c satisfies the GL description of conventional superconductors in the dirty limit l denotes the electron mean free path and in the dirty limit ξ ≈ l . ξ 0 = √ φ 0 /2πH c2, 0 can be calculated from the WHH analysis, ξ 0, x=0.35 = (6.30 ± 0.01) nm and ξ 0, x=0.71 = (6.80 ± 0.01)nm. Fitting Experimental values of ξ(T) and H c1 (T) can be used to determine (T) by using the relation µ 0 H c1 = φ 0 /4π 2 ln( /ξ ) . Knowledge of (T) can provide valuable insights into the nature of superconductivity in the HEA films. We have quantitatively analyzed −2 (T) of both HEA films shown in Fig. 3. To this end, we utilize the BCS superfluid density model in the dirty limit 30 , where k B and �(T) denote Boltzman's constant and the temperature-dependent superconducting quasiparticle gap, respectively.
Our angle-dependent H c1 measurements shown in Fig. 1d reveal an isotropic response of the superconducting state to an external magnetic field rotated in the film plane. This observation supports an isotropic www.nature.com/scientificreports/ superconducting order parameter symmetry, because anisotropic order parameters, such as p-and d-wave, would result in an angle-dependent diamagnetic response. Therefore, we assume an s-wave superconducting order parameter and single-band pairing for our analysis using Eq. (1). The corresponding interpolating BCS gap function reads �(T) = � 0 tanh(α √ T C /T − 1) with 0 = 1.764k B T C . Using this model, we can accurately fit −2 (T) at both alloy compositions as shown in Fig. 3a. Fitting results in zero temperature penetration depths x=0.35 (0) = (896 ± 4) nm and x=0.71 (0) = (245 ± 2) nm. The calculated quasiparticle gaps, which were used for fitting −2 (T) , are displayed in the corresponding insets. The −2 (T) of both HEA films with different compositions is in agreement with the weak-coupling BCS limit, α x=0.35 = α x=0.71 = 1.74 . It is worth noting that we also observe relatively good agreement between experiment and model at intermediate coupling α = 2.2 (dashed lines in Fig. 3), which is consistent with previous reports from heat capacity measurements 20 . Overall, our analysis shows that the superconducting state of (TaNb) 1−x (ZrHfTi) x HEAs can be described with the BCS model for single-band isotropic s-wave superconductivity.
The large degree of compositional disorder is expected to result in a significant on-site potential disorder at the atomic scale 31 . Therefore, the strong dependence of T c on x, T c, x=0.35 = (6.7 ± 1.1) K and T c, x=0.71 = (4.3 ± 1.1) K, invites speculation on the role of disorder for the T c amplitude [21][22][23] . However, our H c2 (T) analysis shown in Fig. 2b reveals a comparable mean free path on the order of 5 to 6 nm at both alloy compositions (see Fig. 2b). This observation is consistent with an almost equivalent mixing entropy. It follows that the microscopic disorder is expected to be of comparable strength in both films, despite their rather different elemental composition. Hence, our measurement results cannot support a disorder-driven mechanism as the origin of the observed T c variations. It is more likely that the T c variations arise from changes to the density of states at the Fermi level induced by electronic doping within a classical BCS framework as reported previously 20 . www.nature.com/scientificreports/

Conclusion
We have experimentally studied the superconducting state of the HEA (TaNb) 1−x (ZrHfTi) x In thin film form with x = 0.35 and x = 0.71 by measuring H c1 (T) and H c2 (T) . Our analysis of (T) is in quantitative agreement with the BCS theory of an isotropic single band s-wave superconductor in the weak coupling limit. The analysis of ξ(T) reveals a comparable amount of disorder at both compositions, l x=0.35 ≈ l x=0.71 . Therefore, we can exclude that the observed variations in T c originate from a disorder-driven mechanism. Further theoretical and experimental studies will be needed to characterize the low-energy electronic structure at various alloy compositions and its influence on T c . Looking ahead, results of such efforts may inform pathways for realizing HEAs with enhanced superconducting T c . Employing penetration-depth measurements to other superconducting HEAs 8 , such as those crystallizing on the CsCl-type lattice, it will be interesting to test whether weak coupling s-wave superconductivity is a common occurrence in these material systems. While the T c variations of bulk superconductivity appear to be independent from disorder, the study of these or other HEAs films in the two-dimensional limit with maximized on-site disorder could offer avenues for exploring T c enhancements through multifractal eigenstates near a quantum critical point 22,23,32 .